Shear stresses of colloidal dispersions at the glass transition in equilibrium and in flow 
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We consider a model dense colloidal dispersion at the glass transition, and investigate the con- 
nection between equilibrium stress fluctuations, seen in linear shear moduli, and the shear stresses 
under strong flow conditions far from equilibrium, viz. flow curves for finite shear rates. To this 
purpose thermosensitive core-shell particles consisting of a polystyrene core and a crosslinked poly(N- 
isopropylacrylamide)(PNIPAM) shell were synthesized. Data over an extended range in shear rates 
and frequencies are compared to theoretical results from integrations through transients and mode 
coupling approaches. The connection between non-linear rheology and glass transition is clarified. 
While the theoretical models semi-quantitatively fit the data taken in fiuid states and the pre- 
dominant elastic response of glass, a yet unaccounted dissipative mechanism is identified in glassy 
states. 
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I. INTRODUCTION 

Complex fluids and soft materials in general are char- 
acterized by a strong variability in their rheological and 
elastic properties under flow and deformationsii Within 
the linear response framework, storage- and loss- (shear) 
moduli describe elastic contributions in solids and dis- 
sipative processes in fluids. Both moduli are connected 
via Kramers-Kronig relations and result from Fourier- 
transformations of a single time-dependent function, the 
shear modulus g^'^{t). Importantly, the linear response 
modulus g'' (t) itself is defined in the quiescent system 
and (only) describes the small shear-stress fluctuations 
always present in thermal equilibriumfii^ 

Viscoelastic materials exhibit both, elastic and dissipa- 
tive, phenomena depending on external control parame- 
ters like temperature and/ or density. The origins of 
the change between fluid and solid like behavior can be 
manifold, including phase transitions of various kinds. 
One mechanism existent quite universally in dense sys- 
tems is the glass transition, that structural rearrange- 
ments of particles become progressively slower It is ac- 
companied by a structural relaxation time which grows 
dramatically. Maxwell was the first to describe this fluid- 
solid transition phenomcnologically. Dispersions consist- 
ing of colloidal, slightly polydisperse (near) hard spheres 
arguably constitute one of the most simple viscoelastic 
systems, where a glass transition has been identified. 
It has been studied in detail by dynamic light scatter- 
ing measurements 16,7,8 ,9 , lo, 1 1 pgnfocal microscopy;^ 
linear^ and non-linear rheologySiii^^^i^i^ 
Computer simulations are available also i^^i^^'^^ Mode 
coupling theory (MCT) has provided a semi-quantitative 
explanation of the observed glass transition phenom- 



ena, albeit neglecting ageing effectssS and decay pro- 
cesses at ultra-long times that may cause (any) colloidal 
glass to flow ultimately4i21i^ Importantly, MCT predicts 
a purely kinetic glass transition and describes it using 
only equilibrium structural input, namely the equilib- 
rium structure factor S^^^ measuring thermal density 
fluctuations. 

The stationary, nonlinear rheological behavior un- 
der steady shearing provides additional insight into the 
physics of dense colloidal dispersions A priori it is 
not clear, whether the mechanisms relevant during glass 
formation also dominate the nonlinear rheology. Sol- 
vent mediated interactions (hydrodynamic interactions), 
which do not affect the equilibrium phase diagram, may 
become crucially important. Also, shear may cause or- 
dering or layering of the particles 1^ Simple phenomcno- 
logical relations between the frequency dependence of the 
linear response and the shear rate dependence of the non- 
linear response, like the Cox-Merz rule, have been for- 
mulated, but often lack firm theoretical support or are 
limited to special shear historiesFii^ 

On the other hand, within a number of theoretical 
approaches a connection between steady state rheology 
and the glass transition has been suggested. Brady 
worked out a scaling description of the rheology based 
on the concept that the structural relaxation arrests 
at random close packing 1^ In the soft glassy rheology 
model, the trap model of glassy relaxation by Bouchaud 
was generalized to describe mechanical deformations and 
ageing i^^i^i^^ The mean field approach to spin glasses 
was generalized to systems with broken detailed balance 
in order to model flow curves of glasses under shear.— i^i 
The application of these novel approaches to colloidal 
dispersions has lead to numerous insights, but has been 
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hindered by the use of unknown parameters in the ap- 
proaches. MCT, also, was generalized to include effects 
of shear,— i^i^ and, within the integrations through tran- 
sients (ITT) approach, to quantitatively describe all as- 
pects of stationary states under steady shearing 
Some aspects of the ITT approach to flow curves have 
already been tested,— but the connection, central in 
the approach, between fluctuations around equilibrium 
and the nonlinear response, has not been investigated 
experimentally up to now. 

In the present contribution we explore the connec- 
tion between structural relaxation close to glassy ar- 
rest and the rheological properties far from equilibrium. 
Thereby we crucially test the ITT approach, which aims 
to unify the understanding of these phenomena. It re- 
quires, as sole input, information on the equilibrium 
structure (namely 5*^), and, first gives a formally ex- 
act generalization of the shear modulus to finite shear 
rates, g{t^j), which is then approximated in a consis- 
tent way. We investigate a model dense colloidal dis- 
persion at the glass transition, and determine its linear 
and nonlinear rheology. Thcrmoscnsitive core-shell par- 
ticles consisting of a polystyrene core and a crosslinked 
poly(N-isopropylacrylamide)(PNIPAM) shell were syn- 
thesized and their dispersions characterized in detail.— 
Data over an extended range in shear rates and frequen- 
cies are compared to theoretical results from MCT and 
ITT. 

The paper is organized as follows: section II summa- 
rizes the equations of the microscopic ITT approach in 
order to provide a self-contained presentation of the the- 
oretical framework. In section III some of the univer- 
sal predictions of ITT are discussed in order to describe 
the phenomenological properties of the non-equilibrium 
transition studied in this work. Building on the univer- 
sal properties, section IV introduces a simplified model 
which reproduces the phenomenology. Section V intro- 
duces the experimental system. Section VI contains the 
main part of the present work, the comparison of com- 
bined measurements of the linear and non-linear rheology 
of the model dispersion with calculations in microscopic 
and simplified theoretical models. A short summary con- 
cludes the paper, while the appendix presents an exten- 
sion of the simplified model used in the main text. 



II. MICROSCOPIC APPROACH 

We consider N spherical particles with radius Rh dis- 
persed in a volume V of solvent (viscosity 77^) with im- 
posed homogeneous, and constant linear shear-flow. The 
flow velocity points along the x-axis and its gradient 
along the y-axis. The motion of the particles (with po- 
sitions ri{t) for 1 = 1,..., N) is described by N coupled 
Langevin cquations^S 



Solvent friction is measured by the Stokes friction coeffi- 
cient C = ^ttt^sRh- The N vectors = —d/dYiU{{rj}) 
denote the interparticle force on particle i deriving from 
potential interactions with all other particles; U is the po- 
tential energy which depends on all particles' positions. 
The solvent shear-flow is given by v^°'^(r) = jyx, and 
the Gaussian white noise force satisfies (with a, P denot- 
ing directions) 

{ft it) //(<')> = ^CksTSo^p % S{t - t') , 

where fc^T is the thermal energy. Each particle expe- 
riences interparticle forces, solvent friction, and random 
kicks. Interaction and friction forces on each particle bal- 
ance on average, so that the particles are at rest in the 
solvent on average. The Stokesian friction is proportional 
to the particle's motion relative to the solvent flow at its 
position; the latter varies linearly with y. The random 
force on the level of each particle satisfies the fluctuation 
dissipation relation. 

An important approximation in Eq. ([T]) is the ne- 
glect of hydrodynamic interactions, which would arise 
from the proper treatment of the solvent flow around 
moving particles i^^^^ In the following we will argue that 
such effects can be neglected at high densities where in- 
terparticle forces hinder and/or prevent structural rear- 
rangements, and where the system is close to arrest into 
an amorphous, metastable solid. Another important ap- 
proximation in Eq. ([T]) is the assumption of a given, 
constant shear rate 7, which does not vary throughout 
the (infinite) system. We start with this assumption 
in the philosophy that, first, homogeneous states should 
be considered, before heterogeneities and confinement ef- 
fects arc taken into account. All difficulties in Eq. ([T|) 
thus arc connected to the many-body interactions given 
by the forces F^, which couple the TV Langevin equa- 
tions. In the absence of interactions, Fj = 0, Eq. ([T|) 
leads to super-diffusive particle motion termed 'Taylor 
dispersion'.— 

While formulation of the considered microscopic model 
handily uses Langevin equations, theoretical analysis 
proceeds more easily from the reformulation of Eq. ([T|) as 
Smoluchowski equation. It describes the temporal evolu- 
tion of the distribution hmction ^'({ri} , t) of the particle 
positions 

dt^{{vi},t) = n^{{v,},t) , (2) 



employing the Smoluchowski operator 
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C(^-v^°'^(r,))=F, + f,: 



(1) 



built with the (bare) diffusion coefficient Dq = ksT/C, 
of a single particle. We assume that the system re- 
laxes into a unique stationary state at long times, so that 
\E'(t — > 00) = ^ s holds. Homogeneous, amorphous sys- 
tems are studied so that the stationary distribution func- 
tion is translationally invariant but anisotropic. Ne- 
glecting ageing, the formal solution of the Smoluchowski 
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FIG. 1: Shear advection of a fluctuation with initial wavevec- 
tor in a;-direction, q{t = 0) = q (1,0,0)^ , and advected 
wavevector at later time q(t>0) — q (1, — 0)"^. At all times, 
q(t) is perpendicular to the planes of constant fluctuation 
amplitude. Note that the magnitude q{t) = q^/i + (7^)^ in- 
creases with time. Brownian motion, neglected in this sketch, 
would smear out the fluctuation. 

equation within the ITT approach can be brought into 
the formiiii^ 

• poo 

= *e + 7^ dt^e l^y c"'* , (4) 
I^B-L Jo 

where the adjoint Smohichowski fjt operator arises from 
partial integrations. It acts on the quantities to be av- 
eraged with '^s- denotes the equilibrum canonical 
distribution function, oc e~^/'^^^'^\ which is the time- 
independent solution of Eq. ([2]) for 7 = 0; in Eq. ([4]), 
it gives the initial distribution at the start of shearing 
(at i = 0). The potential part of the stress tensor 

<7xy = ~Y,i=iFi Vi entered via Vt^e = '^(Txy'^e- The 
simple, exact result Eq. (jl]) is central to the ITT ap- 
proach as it connects steady state properties to time in- 
tegrals formed with the shear-dependent dynamics. The 
latter contains slow intrinsic particle motion. 

In ITT, the evolution towards the stationary distribu- 
tion at infinite times is approximated by following the 
slow structural rearrangements, encoded in the transient 
density correlator $q(i). It is defined by^i^ 

It describes the fate of an equilibrium density fluctua- 
tion with wavevector q, where = X^jLi e^'^ ''^, under 
the combined effect of internal forces, Brownian motion 
and shearing. Note that because of the appearance of 
in Eq. ([U, the average in Eq. (O can be evaluated with 
the equilibrium canonical distribution function, while the 
dynamical evolution contains Brownian motion and shear 
advection. The normalization is given by Sq the equilib- 
rium structure factori2i2^ for wavevector modulus g = |q|. 
The advected wavevector enters in Eq. ^ 

q(0 = q - 7* 9^ y , (6) 

where unit-vector y points in y-direction) The time- 
dependence in q(t) results from the affine particle mo- 



tion with the shear flow of the solvent. Translational 
invariance under shear dictates that at a time t later, 
the equilibrium density fluctuation Jp* has a nonvanish- 
ing overlap only with the advected fluctuation (5pq(t); see 
Fig. [U where a non-decorrelating fluctuation is sketched 
under shear. In the case of vanishing Brownian motion, 
viz. £)o = in Eq. ([3]), we find $q(t) = 1, because 
the advected wavevector takes account of simple affine 
particle motion4^ The relaxation of <&q(t) thus heralds 
decay of structural correlations. Within ITT, the time 
integral over such structural decorrelations provides an 
approximation to the stationary state: 

*s « *e+ (7) 

with S'f. = dSk/dk^ The last term in brackets in Eq. 
([7|) expresses, that the expectation value of a general 
fluctuation A in ITT-approximation contains the (equi- 
librium) overlap with the local structure, {g^^ gk A)'''^^^\ 
The difference between the equilibrium and stationary 
distribution functions then follows from integrating over 
time the spatially resolved (viz. wavevector dependent) 
density variations. 

The general results for s , the exact one of Eq. (|4|) and 
the approximation Eq. ([7]), can be applied to compute 
stationary expectation values like for example the ther- 
modynamic transverse stress, cr(7) = {axy)/V. Equation 
(HI) leads to an exact non-linear Green-Kubo relation: 

/•OO 

a{^)^j dtg{t,i), (8) 
Jo 

where the generalized shear modulus git^j) depends on 
shear rate via the Smoluchowski operator from Eq. ([3]) 

5(^,7) = ^^(^..e"'*^..>^^=°^ (9) 

In ITT, the slow stress fluctuations in .9(^,7) are ap- 
proximated by following the slow structural rearrange- 
ments, encoded in the transient density correlators. The 
generalized modulus becomes, using the approximation 
Eq. (O, or, equivalently, performing a mode coupling 
approximation)^2illd^ 

Summation over wavevcctors has been turned into inte- 
gration in Eq. (jlO|) considering an infinite system. 

The familiar shear modulus of linear response the- 
ory describes thermodynamic stress fluctuations in equi- 
librium, and is obtained from Eqs. (|9ll0p by setting 
7 = 0.-'^'^ While Eq. ^ then gives the exact Green- 
Kubo relation, the approximation Eq. pH)) turns into the 
well-studied MCT formula. For finite shear rates, Eq. 
(fTU)) describes how affine particle motion causes stress 
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fluctuations to explore shorter and shorter length scales. 
There the effective forces, as measured by the gradient of 
the direct correlation function, S'j^/Sl = ncj,. = ndck/dk, 
become smaller, and vanish asympotically, c'i^_^^ — > 0; 
the direct correlation function Ck is connected to the 
structure factor via the Ornstein-Zernicke equation Sk = 
1/(1 — ncfe), where n = A^/l^isthe particle density. Note , 
that the equilibrium structure suffices to quantify the ef- 
fective interactions, while shear just pushes the fluctua- 
tions around on the 'equilibrium energy landscape'. 

Structural rearrangements of the dispersion affected 
by Brownian motion is encoded in the transient density 
correlator. Shear induced affine motion, viz. the case 
Do = 0, is not sufficient to cause $k(i) to decay. Brown- 
ian motion of the quiescent correlator ^l^~^\t) leads at 
high densities to a slow structural process which arrests 
at long times in (metastable) glass states. Thus the com- 
bination of structural relaxation and shear is interesting. 
The interplay between intrinsic structural motion and 
shearing in <I>k(0 captured by (i) first a formally exact 
Zwanzig-Mori type equation of motion, and (ii) second a 
mode coupling factorisation in the memory function built 
with longitudinal stress fluctuations The equation of 
motion for the transient density correlators is 

at*q(i)+rq(t) |*qW+_^ di'mq(t,t')5t'$q(t')| =0, 

(11) 

where the initial decay rate Tq{t) = Doq^{t)/ Sq^t) gen- 
eralizes the familiar result from linear response theory to 
advected wavevectors; it contains Taylor dispersion. The 
memory equation contains fluctuating stresses and simi- 
larly like g{t, 7) in Eq. ([7]), is calculated in mode coupfing 
approximation 

™q(i. = ^ E ^qkp(i, t') <I>k(t') {t - t') <i>p(,,) {t - 1') , 
k 

(12) 

where we abbreviated p = q — k. The vertex generalizes 
the expression in the quiescent casei^ 

V (+ i'\ -^qlt) '^k(f') 'S'p(f/) 

Kqkp(t,t ) = g2(-^) g2(^/-) l^qkp(l) Vqkp(t ) , 

Vqkp(t) = q(t) • (k(t) riCi,(t) +p(t) ncp(t)) . (13) 

With shear, wavevectors in Eq. are advected accord- 
ing to Eq. dSl) 

Equations (|7|llll2p . with the specific example of the 
generalized shear modulus Eq. PU)) . form a closed set of 
equations determining rheological properties of a sheared 
dispersion from equilibrium structural input Only 
the static structure factor Sq is required to predict (i) 
the time dependent shear modulus within linear response, 
^''(i) — di^ii ~ 0)j ^^^id (m) the stationary stress (t(7) 
from Eq. ([8]). The loss and storage moduli of small 
amplitude oscillatory shear measurcmcntsiiiSi follow from 
Eq. ^ in the linear response case {i) 

/>oo 

G'[uj)+iG"{uo)^iu / (ite-*"*g(i,7 = 0) . (14) 
Jo 



While, in the linear response regime, modulus and den- 
sity correlator are measurable quantities, outside the lin- 
ear regime, both quantities serve as tools in the ITT ap- 
proach only. The transient correlator and shear modu- 
lus provide a route to the stationary averages, because 
they describe the decay of equilibrium fluctuations un- 
der external shear, and their time integral provides an 
approximation for the stationary distribution function, 
see Eq. Determination of the frequency dependent 
moduli under large amplitude oscillatory shear has be- 
come possible recently only;^ and requires an extension 
of the present approach to time dependent shear rates in 
Eq. 



III. UNIVERSAL ASPECTS 

The summarized microscopic ITT equations contain 
a bifurcation in the long-time behavior of ^qif), which 
corresponds to a non-equilibrium transition between a 
fluid and a shear-molten glassy state; it is described in 
this section. Close to the transition, (rather) universal 
predictions can be made about the non-linear dispersion 
rhcology and the steady state properties. The central 
predictions are introduced in this section and summa- 
rized in the overview figure [21 It is obtained from the 
schematic model which is also used to analyse the data, 
and which is introduced in the following Sect. IV. 

A dimensionless separation parameter e measures the 
distance to the transition which is situated at e = 0. 
A fluid state (e < 0) possesses a (Newtonian) vis- 
cosity, ?7o(£ < 0) = lim.y^o o'(7)/7, and shows shear- 
thinning upon increasing 7. Via the relation ryo = 
lim;^_»o G"{u})/uj, the Newtonian viscosity can also be 
taken from the loss modulus at low frequencies, where 
G"{uj) dominates over the storage modulus. The latter 
varies like G'{u! ^ 0) ^ uj'^. A glass (e > 0), in the 
absence of flow, possesses an elastic constant Goo , which 
can be measured in the elastic shear modulus G'{lu) in 
the limit of low frequencies, G"(ci; 0, e > 0) — ^ Goois)- 
Here the storage modulus dominates over the loss one, 
which drops like G"{uj — > 0) ~ w. (Note that the high 
frequency modulus G'^ = G"((jJ ^ 00) is characteristic of 
the particle interactions and exists in fluid and solid 
states.) Enforcing steady shear flow melts the glass. The 
stationary stress of the shear-molten glass always exceeds 
a (dynamic) yield stress. For decreasing shear rate, the 
viscosity increases like I/7, and the stress levels off onto 
the yield-stress plateau, (7(7 ^ 0,e > 0) ^ '^'^{s)- 

Close to the transition, the zero-shear viscosity 779, 
the clastic constant Goo, and the yield stress show 
universal anomalies as functions of the distance to the 
transition: the viscosity diverges in a power-law rjQ^e — *■ 
0— ) ~ (— e)^''' with material dependent exponent 7 
around 2 — 3, the elastic constant increases like a square- 
root Goo(£ — > 0-I-) — G^ ~ y^, and the dynamic yield 
stress ^^(e 0+) also increases with infinite slope above 
its value at the bifurcation. The quantities G^ and 
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FIG. 2: Overview of the properties of the F^j niodel char- 
acteristic for the transition between fluid and yielding glass. 
The upper panel shows numerically obtained transient cor- 
relators for e = 0.01 (black curves), £ = (red), 
e = -0.005 (green), and e = -0.01 (blue). The shear 
rates are |7/r| = (thick solid lines), |7/r| = lO"** (dot- 
ted lines), and |7/r| — 10~^ (dashed lines). For the glass 
state at £ = 0.01 (black), |7/r| = 10"* (dashed-dotted-line) 
is also included. All curves were calculated with 7c = 0.1 
and ?7oo = 0. The thin solid lines give the factorization re- 
sult Eq. I|15|) with scaling functions Q for |7/r| = 10"®; la- 
bel a marks the critical law (|17|l . and label b marks the von 
Schweidler-law (|18p . The critical glass form factor fc is indi- 
cated. The inset shows the flow curves for the same values 
for e. The thin black bar shows the yield stress cr^ for e = 0. 
The lower panel shows the viscoelastic storage (solid line) and 
loss (broken line modulus for the same values of e. The thin 
green lines are the Fourier-transformed factorization result 
Eq. (|15|l with scaling function Q taken from the upper panel 
for e = —0.005. The dashed-dotted lines show the fit for- 
mula Eq. (|29|l for the spectrum in the minimum-region with 
Gminlv„ = 0.0262, t^n,in/r = 0.000457 at e = -0.005 (green) 
and Gmin/va = 0.0370, t^min/F = 0.00105 at e = -0.01 
(blue). The elastic constant at the transition is marked 
also, while the high frequency asymptote G'^x, ~ G'iui — > oo) 
is not labeled explicitly. 



(T^ denote the respective values at the transition point 
e = 0, and measure the jump in the elastic constant and 
in the yield stress at the glass transition; in the fluid 
state, Gc»(£ < 0) = and (T+(e < 0) = hold. 

The described results follow from the stability analy- 
sis of Eqs. (|llll2p around an arrested, glassy structure 



fq of the transient correlator 4^*^ Considering the time 
window where *&q(i) is metastable and close to arrest at 
/q, and taking all control parameters like density, tem- 
perature, etc. to be close to the values at the transition, 
the stability analysis yields the 'factorization' between 
spatial and temporal dependences 



= fa 



hq g{t/tQ,s,'yto ) 



(15) 



where the (isotropic) glass form factor and criti- 
cal amplitude hq describe the spatial properties of the 
metastable glassy state. The critical glass form factor /| 
gives the long-lived component of density fluctuations, 
and hq captures local particle rearrangements. Both can 
be taken as constants independent on shear rate and den- 
sity, as they are evaluated from the vertices in Eq. (fT3l) 
at the transition point. All time-dependence and (sensi- 
tive) dependence on the external control parameters is 
contained in the function Q, which often is called '/?- 
correlator' and obeys the non-linear stability equation 



^S\t) = ^ fdt' g{t-t')g{t') , (16) 

at Jo 



with initial condition 



Git ^ 0) ^ {t/t„y 



(17) 



The two parameters A and c'^-' in Eq. ([TB|) are deter- 
mined by the static structure factor at the transition 
point, and take values around A w 0.73 and c^^' ~ 3 
for Sq taken from Percus-Yevick approximation^ for hard 
sphere interactions ^ii^i^ The transition point then lies 
at packing fraction (pc = ^n-cR'^j ~ 0.52 (index c for 
critical) , and the separation parameter measures the rel- 
ative distance, e = C {(/) — (t)c)l4>c with C « 1.3. The 
'critical' exponent a is given by the exponent parameter 
A via A = r(l - a)Vr(l - 2a)^ 



The time scale in Eq. (|T7|) provides the means to 
match the function Q{t) to the microscopic, short-time 
dynamics. The Eqs. pill2p contain a simplified descrip- 
tion of the short time dynamics in colloidal dispersions 
via the initial decay rate rq(i). From this model for the 
short-time dynamics, the time scale ig ~ 1.6 10"^i?|^/£'o 
is obtained. Solvent mediated effects on the short time 
dynamics are well known and are neglected in Tq(t) in 
Eq. pT|) . Within the ITT approach, one finds that if hy- 
drodynamic interactions were included in Eq. (jlip , all of 
the mentioned universal predictions would remain true. 
Only the value of to will be shifted and depend on the 
short time hydrodynamic interactions. For the quiescent 
glass transition this has been discussed within MCT^ 
and ITT extends this to driven cases. This statement 
remains valid in ITT, as long as the hydrodynamic in- 
teractions do not affect the mode coupling vertex in Eq. 
(fT^ . In this sense, hydrodynamic interactions can be 
incorporated into the theory of the glass transition, and 
amount to a rescaling of the matching time toi only. 

Obviously, the matching time to a-lso provides an upper 
cut-off for the time window of the structural relaxation. 
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At times shorter than to the specific short-time dynamics 
matters. The condition 'ytQ ^ 1 foUows and translates 
into a restriction for the accessible range of shear rates, 
7^7*, where the upper-cut off shear rate 7, is con- 
nected to the matching time. 

The parameters e, A and c'^^ in Eq. p6|) can be deter- 
mined from the equilibrium structure factor Sq at or close 
to the transition, and, together with and the shear rate 
7 they capture the essence of the rheological anomalies in 
dense dispersions. A divergent viscosity follows from the 
prediction of a strongly increasing final relaxation time 
in in the quiescent fluid phase 

g{t ^oo,e<0,j^O)~^~ {t/Tf , with — cx {^ey . 

(18) 

The entailed temporal power law, termed von Schweidler 
law, initiates the final decay of the correlators, which 
has a density and temperature independent shape ^q{i). 
The final decay, often termed a-relaxation, depends on 
s only via the time scale T{e) which rescales the time, 
t = t/r. Equation (fTB|) establishes the crucial time scale 
separation between to and r, the divergence of t, and 
the stretching (non-exponentiality) of the final decay; it 
also gives the values of the exponents via A = r(l + 
6)Vr(l-h2&), and 7 = {l/a + l/b)/2. Using Eq. the 
divergence of the Newtonian viscosity foUows^i^ During 
the final decay the shear modulus becomes a function of 
rescaled time, g{i = t/T,'j = 0), leading to t]o oc T{e); 
its initial value is given by the elastic constant at the 
transition, g{i <cl,£^0— ,7 = 0) = G^. 

On the glassy side of the transition, e > 0, the tran- 
sient density fluctuations stays close to a plateau value 
for intermediate times which increases when going deeper 
into the glass, 



g(to<t<l/|7|,e>0)^^Y-^+0(e). (19) 

Entered into Eq. ([TU]) . the square-root dependence of the 
plateau value translates into the square-root anomaly of 
the elastic constant Goo, and causes the increase of the 
yield stress close to the glass transition. 

Only, for vanishing shear rate, 7 = 0, an ideal glass 
state exists in the ITT approach for steady shearing. All 
density correlators arrest at a long time limit, which from 
Eq. (|19p close to the transition is given by $q(t — > 00, e > 
0,7 = 0) = /g = fq+hq^e/{l - X)+0{e). Consequently 
the modulus remains elastic at long times, g{t 00, e > 
0,7 = 0) = Goc > 0. Any (infinitesimal) shear rate, 
however, melts the glass and causes a final decay of the 
transient correlators. The function Q initiates the decay 
around the critical plateau of the transient correlators 
and sets the common time scale for the final decay under 
shear 



g{t^^,e>0)^-Jj-j\jt\. (20) 



Under shear all correlators decay from the plateau as 
function of {jtl. Steady shearing thus prevents non- 
ergodic arrest and restores ergodicity. This aspect of 
Eq. (|16p has two important ramifications for the steady 
state of shear molten glasses4ii^ First, ITT finds that 
shear melts a glass and produces a unique steady state 
at long times. This conclusion is restricted by the as- 
sumption of homogeneous states and excludes the possi- 
ble existence of ordering or layering under shear. Also, 
ageing was neglected, which could remain because of non- 
ergodicity in the initial quiescent state. (Ergodicity of the 
sheared state however suggests ageing to be unimportant 
under shear. ^^i^^ ) Second, all stationary averages, which 
in ITT are obtained from integrating up the transient 
fluctuations, do not exhibit a linear response regime in 
the glass. Rather they take finite values for vanishing 
shear rate, which jump discontinuously at the glass tran- 
sition. This holds because the shear-driven decay of Eq. 
(PP)) initiates a scaling law where the transient correla- 
tors decay as function of |7t| down from the plateau fq 
to zero, denoted as $q(t|7|). When entered into Eq. 
([7]), time appears only in the combination together with 
shear rate and thus after time integration the shear rate 
dependence drops out, yielding a finite result even in the 
limit of infinitesimal shear rate. Prominent example of 
a stationary average that has no linear response regime 
with respect to 7 in the glass phase is the shear stress 
"■(i'lE > 0). It takes finite values for vanishing shear 
rate, (J'^{s) = 17(7 ^ 0,e > 0), and jumps at the glass 
transition from zero to a finite value. Because of Eq. (fTO|) 
it increases rapidly when moving deeper into the glass. 



IV. SCHEMATIC MODEL 

The universal aspects described in the previous sec- 
tion are contained in any ITT model that contains the 
central bifurcation scenario and recovers Eqs. pSIlGp . 
Equation (|15p states that spatial and temporal depen- 
dences decouple in the intermediate time window. Thus 
it is possible to investigate ITT models without proper 
spatial resolution. Because of the technical difficulty to 
evaluate the anisotropic functionals in Eqs. (|10I12|) . it is 
useful to restrict the description to few or to a single tran- 
sient correlator. In the schematic F^2^ -model a single 
'typical' density correlator $(t), conveniently normalized 
according to ^{t 0) = 1 — Tt, obeys a Zwanzig-Mori 
memory equation which is modeled according to Eq. pip 

at$(t)+r|$(t) +^ dt' m{t-t') dt'^it')^ ^0 . (21) 

The parameter F mimics the microscopic dynamics of 
the 'typical' density correlator chosen in Eq. (pij) . and 
will depend on structural and hydrodynamic correla- 
tions. The memory function describes stress fluctuations 
which become more sluggish together with density fluctu- 
ations, because slow structural rearrangements dominate 
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all quantities. A self consistent approximation closing 
the equations of motion is made mimicking Eq. (|12p . In 

the F^^^ -model one includes a linear term (absent in Eq. 
(fT2| ) in order to {i) sweep out the full range of A values 
in Eq. p6p . and (ii) retain algebraic simplicity: 



m{t) 



1 + iit/jcf 



(22) 



This model, for the quiescent case 7 = 0, had been 
suggested by Gotze in 1984f2L^ and describes the de- 
velopment of slow structural relaxation upon increasing 
the coupling vertices Vi > 0; they mimic the dependence 
of the vertices in Eq. p2[) at 7 = on the equilibrium 
structure given by Sq. Under shear an explicit time de- 
pendence of the couplings in m{t) captures the accel- 
erated loss of memory by shear advection in Eq. (H^). 
Shearing causes the dynamics to decay for long times, 
because fluctuations are advected to smaller wavelengths 
where small scale Brownian motion relaxes them. Equa- 
tions (|2TT22)1 lead, with $(t) = + (1 - /=)2 ^(i, e, 7), 
and the choice of the vertices V2 ~ V2 = 2, and vi = 
Vi + e {1 — f^)/ where = 0.828, to the critical glass 
form factor = 0.293 and to the stability equation p^ . 
with parameters 

A = 0.707, c^^^ = 0.586/7^^ , and to = 0.426/r . 

The choice of transition point {vi, V2) is motivated by its 
repeated use in the literature. Actually, there is a line of 
glass transitions where the long time limit / = ^{t 00) 
jumps discontinuously. It is parameterized by (f J, u^) = 
((2A - 1), 1)/A2 with 0.5 < A < 1, and f^l-X. The 
present choice is just a typical one, which corresponds 
to the given typical A-value. The separation parameter 
£ is the crucial control parameter as it takes the system 
through the transition. The parameter 7c is a scale for 
the magnitude of strain that is required in order for the 
accumulated strain -yt to matter.— In Eq. (fT6|) . it is 
connected to the parameter c^'^^ . 

For simplicity, the quadratic dependence of the gener- 
alized shear modulus on density fluctuations is retained 
from the microscopic Eq. ((T0|) . It simplifies because only 
one density mode is considered, and as, for simplicity, a 
possible dependence of the vertex (prefactor) on shear 
is neglected 



g{t)^v,'S>\t) + 7^oo Sit) . 



(23) 



The parameter rjaa characterizes a short-time, high fre- 
quency viscosity and models viscous processes which re- 
quire no structural relaxation. Together with T (respec- 
tively to), it is the only model parameter affected by 
solvent mediated interactions. Steady state shear stress 
under constant shearing, and viscosity then follow via 
integrating up the generalized modulus: 



(7 = 77 7 = 7 



digit) =^ / dt v„^\t) + -j ri^ 



Also, when setting shear rate 7 = in Eqs. (|21l22p . so 
that the schematic correlator belongs to the quiescent, 
equilibrium system, the frequency dependent moduli are 
obtained from Fourier transforming: 



G'iuj)+iG"iu;) ^iio / e" 
Jo 

(25) 

Because of the vanishing of the Fourier-integral in Eq. 
(|25p for high frequencies, the parameter r^oo can be iden- 
tified as high frequency viscosity: 



lim G"(w)/w = ry^ 



with rj'^ 



(26) 



At high shear, on the other hand, Eq. leads to a 

vanishing of 7n(i), and Eq. (j2ip gives an exponential 
decay of the transient correlator, $(t) —* e^'"* for 7^0. 
The high shear viscosity thus becomes 



7— >c» Z i 



2r 



(27) 



Representative solutions of the F^^''-model are sum- 
marized in Figl^l which bring out the discussed universal 
aspects included in all ITT models. 



V. EXPERIMENTAL SYSTEM AND METHODS 

The particles consist of a solid core of poly(styrcne) 
onto which a network of crosslinked poly(N- 
isopropylacrylamide) (PNIPAM) is affixed. The 
degree of crosslinking of the PNIPAM shell effected by 
the crosslinker N,N'-methylenebisacrylamide (BIS) was 
5 Mol %. The core-shell type PS-NIPAM particles were 
synthesized, purified and characterized as described in 
refi^ Immersed in water the shell swells at low tem- 
peratures. Raising the temperature above 32°C leads 
to a volume transition within the shell. To investigate 
the structure and swelling of the particles cryogenic 
transmission electron microscopy and dynamic light 
scattering have been used. 

Screening the remaining electrostatic interactions by 
adding 5.10~^ molL~^ KCl, the system crystallises as 
hard spheres.— Experimental details on the characteri- 
zation of the particles and on the determination of the 
effective volume are given elsewhere4^ The dependence 
of (f>cS on the temperature is given by the hydrodynamic 
radius Rh determined from the dynamic light scattering 
in the dilute regime. Rh was linearly extrapolated be- 
tween 14 and 25°C (i?i/ = -0.85925T + 123.78 with T 
the temperature in °C) as described recently^^ and (j)cs 
was calculated following the relation 



Pes 



Rh 
~R~ 



(28) 



(24) 



with R the radius of the core determined by cryogenic 
transmission electron microscopy (R = 52 nm)^ c the 
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concentration in wt % and k a rescaling constant. In 
order to determine fc, an experimental phase diagram 
has been achieved by determining the crystal fraction of 
the samples from the position of the coexistence liquid- 
crystal boundaries after sedimentation. This was linearly 
extrapolated to identify the beginning and the end of the 
coexistence domain. The experimental phase diagram 
of the suspensions of the core-shell particles was rcscaled 
with the constant k = 0.4814 to the freezing volume frac- 
tion for hard spheres (f>F = 0.494i^ 

Three instruments were employed in the present study 
to investigate the rheological properties of the suspen- 
sions. The flow behaviour and the linear viscoelastic 
properties for the range of the low frequencies were mea- 
sured with a stress-controlled rotational rheometer MCR 
301 (Anton Paar), equipped with a Searle system (cup 
diameter: 28.929 mm, bob diameter: 26.673 mm, bob 
length: 39.997 mm). Measurements have been performed 
on 12 ml solution and the temperature was set with an ac- 
curacy of ± 0.05°C. The shear stress a versus the shear 
rate 7 (flow curve) was measured after a pre-shearing 
of 7 = 100 for two minutes and a timesweep of 1 
hour at 1 Hz and 1 % deformation in the linear regime. 
The flowcurves experiment were performed setting 7, first 
with increasing 7 from 7 = 10~^— lO'^ with a logarith- 
mic time ramp from 600 to 20 s, and then with decreasing 
7. The stationarity has been checked by step-flow experi- 
ments in the glassy state for the highest effective volume 
fraction {(peS = 0.622). The frequency dependence of 
the loss G" and elastic G' moduli has been measured 
for 1 % strain from 15 to 10^"^ Hz with a logarithmic 
time ramp from 20 to 600 s. The dependence upon the 
strain has been checked and confirmes that all the mea- 
surements were performed in the linear regime. The fre- 
quency dependence was tested for two different sample 
histories. The experiments were first performed without 
pre-shearing after the timesweep, before the flowcurves 
experiments, and then after the flowcurves experiments 
10 s after two minutes pre-shearing at 7 = 100 s^^ to 
melt eventual crystallites. We only considered experi- 
ments performed after pro-shearing in the following dis- 
cussion of G' and G" for the lowest frequencies. 

Additional rheological experiments were carried out on 
Piezoelectric vibrator (FAV)^ and cylindrical torsional 
resonato r^^i^^ supplied by the Institut fiir dynamische 
Materialpriifung, Ulm. Germany. The PAV was oper- 
ated from 10 to 3000 Hz. The solution is placed be- 
tween two thick stainless steel plates. The upper one 
remains static whereas the lower is cemented to piezo- 
electric elements. The gap was adjusted with a 100 
ring. One set of piezoelectric elements is driven by an 
ac voltage to induce the squeezing of the material be- 
tween the two plates, whereas the second set gives the 
output voltage. Experimental details concerning this in- 
strument are given elsewhere.— Only the measurements 
in the glassy state have been performed with the PAV as 
the instrument does not allow any pre-shearing. 

The cylindrical torsional resonator used was operated 



at a single frequency (26 kHz). The experimental pro- 
cedure and the evaluation of data have been described 
recently.-^^i^ 

The effect of the shear rate 7 on the particle dynam- 
ics is measured by the Peclet number;^ Peg = jRj^/DQ, 
which compares the rate of shear flow with the time an 
isolated particle requires to diffuse a distance identical 
to its radius. Similarly, frequency will be reported in the 
following rescaled by this diffusion time, uj' = ujRj^/Do. 
The self diffusion coefficient Dq at infinite dilution was 
calculated from the hydrodynamic radius Rh and the vis- 
cosity of the solvent 77^ with the Stokes- Einstein relation 
so that Do = ksT /Q-K-qsRH ■ In dense dispersions, how- 
ever, the structural rearrangements proceed far slower 
than diffusion at infinite dilution, and therefore, very 
small Peclet numbers and rescaled frequencies uj' are of 
interest in the following. Stresses will be measured in 
units of ksT/R^ in the following. 

VI. COMPARISON OF THEORY AND 
EXPERIMENT 

Shear stresses measured in non-linear response of the 
dispersion under strong steady shearing, and frequency 
dependent shear moduli arising from thermal shear stress 
fluctuations in the quiescent dispersion were measured 

and fltted with results from the schematic Fj2''"inodel. 
Some results from the microscopic MCT for the equi- 
librium moduli arc included also; see Sect. VI. C for 
more details.^ In the following discussion, we first start 
with more general observations on typical fluid and glass 
like data, and then proceed to a more detailed analysis. 
Figures [3] and [4] show measurements in fluid states, at 
0off = 0.540 and (j)cff = 0.567, respectively, while Fig. [5] 
was obtained in the glass at 4>cs = 0.627. From the 
fits to all (j>cS, the glass transition value (p^g = 0.58 was 
obtained; which agrees well with the measurements on 
classical hard sphere colloids . 



A. Crystallisation effects 

We start the comparison of experimental and theoret- 
ical results by recalling the interpretation of time in the 
ITT approach. Outside the linear response regime, both 
<i>(t) and g{t) describe the decorrelation of equilibrium, 
fluid-like fluctuations under shear and internal motion. 
Integrating through the transients provides the steady 
state averages, like the stress. While theory finds that 
transient fluctuations always relax under shear, real sys- 
tems may either remain in metastable states if 7 is too 
small to shear melt them, or undergo transitions to het- 
erogeneous states for some parameters. In these circum- 
stances, the theory can not be applied, and the rheo- 
logical response of the system, presumably, is dominated 
by the heterogeneities. Thus, care needs to be taken 
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FIG. 3: The reduced flow curves and the corresponding mod- 
uh for a fluid state at 13.01wt%, T= 20°C, and (j)cB = 0.540. 
Flow curves measured proceeding from higher to lower shear 
rates (called 'down' flow curves) and dynamic experiments 
were fitted where effects from crystallisation can be neglected; 
the lower limits of the unafi'ected-data regions are marked by 
vertical bars. The red lines show the fits with the schematic 
'-model while the blue lines show the results from micro- 
scopic MCT (solid G', broken G"), with parameters: e = 



-0.05, 



= 0.15, and rjao ~ Q.^ksT / [DoRh); the moduli 



were scaled up by a factor 



1.4. 



critical shear rate 7cr « 4s~"'^ no hysteresis has been ob- 
served, which proved that all the crystallites have been 
molten. In the present work we focus on the 'down' flow 
curves, and consider only data taken either for 7 > 7cr, 
or (for 7 < 7cr) before the time crystallisation sets in. 
This time was estimated from timesweep experiments as 
the time where crystallisation caused a 10% deviation of 
the complex modulus G* . Thus, only the portions of the 
flow curves unaffected by crystallisation are taken into 
account; in Figs. [3] and [J] vertical bars denote the lim- 
its. We find that the effect of crystallisation on the flow 
curves is maximal around = 0.55 and becomes progres- 
sively smaller and shifts to lower shear rates for higher 
densities; see Figs. [3l|4] and [5l This agrees with the notion 
that the glass transition slows down the kinetics of crys- 
tallisation and causes the average size of crystallites to 
shrink.— For the highest densities, which are in the glass 
without shear, the hysteresis at the lowest 7 can be at- 
tributed to a non stationarity of the up curve; see Fig. [51 
This effect has been confirmed by step flow experiments, 
but does not affect the back curves. 

The linear response moduli similarly are affected by the 
presence of small crystallites at low frequencies. G'{iu) 
and G"{uj) increase above the behavior expected for a so- 
lution {G'{iu — > 0) — ^ ijaU! and G"{u! ^ 0) ^ cuj'^) even 
at low density, and exhibit clastic contributions at low 
frequencies (apparent from G'{uj) > G"{uj)); see Figs. [3] 
andSl This effect follows the crystallisation of the system 
during the measurement after the shearing at 7 = 100 
s^^ . Only data will be considered in the following which 
were collected before the crystallisation time. For higher 
effective volume fraction other effects such as ageing and 
an ultra-slow process had to be taken into account and 
will be discussed more in detail in the next section. 



Shapes of flow curves and moduli and their 
relations 



in experiments, in order to prevent phase transitions, 
and to shear melt arrested structures, before data can 
be recorded. 

The small size polydispersity of the present particles 
enables the system to grow crystallites according to its 
equilibrium phase diagram. Fortunately, when record- 
ing flow curves, viz. stress as function of shear rate, 
data can be taken when decreasing the shear rate. We 
find that the resulting 'down' flow curves correspond to 
amorphous states and reach the expected low-7 asymp- 
totes (cr = rjQ 7, sec Fig. except for very low 7, when 
an increase in stress indicates the formation of crystal- 
lites. 'Up' flow curves, however, obtained when moving 
upwards in shear rate during the measurement of the 
stress are affected by crystallites formed after the initial 
shearing at 100 during the timesweep and the first 
frequency sweep experiments. See the hysteresis between 
'up' and 'down' flow curves in Figs. [3] and [51 where mea- 
surements for two fluid densities arc reported. Above a 



The flow curves and moduli exhibit a qualitative 
change when increasing the effective packing fraction 
from around 50% to above 60%. For lower densities (sec 
Fig. [3]), the flow curves exhibit a Newtonian viscosity t]o 
for small shear rates, followed by a sublinear increase of 
the stress with 7; viz. a region of shear thinning behavior. 
For the same densities, the frequency dependent spec- 
tra exhibit a broad peak or shoulder, which corresponds 
to the final or a-relaxation discussed in Sect. III. Its 
peak position (or alternatively the crossing of the mod- 
uli, G" = G") is roughly given by cot = 1 (see Fig. [J). 
These properties characterize a viscoelastic fluid. For 
higher density, see Fig. [51 the stress in the flow curve 
remains above a finite yield value even for the small- 
est shear rates investigated. The corresponding storage 
modulus exhibits an elastic plateau at low frequencies. 
The loss modulus drops far below the elastic one. These 
observations characterize a soft solid. The loss modulus 
rises again at very low frequencies, which may indicate 
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FIG. 4: The reduced flow curves and the corresponding mod- 
uli (like in Fig. E]) for a fluid state at 13.01wt%, T= 18°C, and 
4'cS ~ 0.567. The vertical bars mark the minimal Peclet num- 
ber or rescaled frequency for which the influence of crystalli- 
sation can be neglected. Microscopic parameters: e — —0.01, 
jf- = 0.15, and 7700 = 0.3 ksT / (DoRh); moduli scale factor 
c/= 1.4. 



that the colloidal solid at this density is metastable and 
may have a finite lifetime (an ultra-slow process is dis- 
cussed in Sect. VI. E). 

Simple relations, like the 'Cox-Merz rule', have some- 
times been used in the past to compare the shapes of the 
flow curves 17(7) with the shapes of the dissipative mod- 
ulus G"{ll!). Both quantities can be interpreted in terms 
of a (generalized) viscosity, on the one hand as function 
of shear rate 77(7) = cr('j)/'y, and on the other hand as 
hmction of frequency ri{Lu) = G"{uj)/uj. The Cox-Merz 
rule states that the functional forms of both viscosities 
coincide. 

Figures [3] to [5] provide a sensitive test of relations in 
the shapes of a{j) and G"{lo). Figure [4] shows most con- 
clusively, that no simple relation between the far-from 
equilibrium stress as function of external rate of shear- 
ing exists with the equilibrium stress fluctuations at the 
corresponding frequency. While 17(7) increases monoton- 
ically, the dissipative modulus G" {lo) exhibits a minimum 
for fluid states close to the glass transition. It separates 
the low-lying final relaxation process in the fluid from the 
higher-frequency relaxation. 

As shown in Fig. [21 the frequency dependence of G" in 



FIG. 5: The reduced How curves and the corresponding mod- 
uli for a glass state at 13.01wt%, T= 14°C, and 0cff = 0.622. 
See figure [3] for further explanations. Microscopic parame- 
ters: e = 0.03, §^ = 0.08, and r?oc = O.ZksT / {DoRh)\ mod- 
uli scale factor Cy = 1.4 (blue). Curves from the schematic 
Fj^^-model with an additional dissipative process included 
(Eg. IA1|) are shown as dashed lines; 5 = lO^'^ F (long dashes, 
light green) and 5 = 10~* F (short dashes, dark green). Here 
F = 88Do/Rji. The red curves give the schematic model 
calculations for identical parameters but without additional 
dissipative process (viz. 5 — 0). 



the minimum region is given by the scaling function Q of 
Sect. Ill, which describes the minimum as crossover be- 
tween two power laws. The approximation for the mod- 
ulus around the minimum 



G"iiu) 



V uj 



(29) 



has been found in the quiescent fluid (e < 0, 7 = 0), 
and works quantitatively if the relaxation time r is large, 
viz. time scale separation holds for small \e\^ The pa- 
rameters in this approximation follow from Eqs. ()17I18|) 
which give Gmin oc and tj,„in oc (— e)^/^°. For pack- 
ing fractions too far below the glass transition, the final 
relaxation process is not clearly separated from the high 
frequency relaxation. This holds in Fig.[3l where the final 
structural decay process only forms a shoulder. Closer to 
the transition, in Fig. [4l it is separated, but crystallisa- 
tion effects prevent us from fitting Eq. ([29]) to the data. 

Asymptotic power-law expansions of a{'y) exist close 
to the glass transition, which were deduced from the sta- 



bility analysis in Sect, ni j i^^i^^ yet we refrain from en- 
tering their detailed discussion and describe the qualita- 
tive behavior in the following. For the same parameters 
in the fluid, where the minimum in G"{lj) appears, the 
flow curves follow a S-shapc in a double logarithmic plot, 
crossing over from a linear behavior cr = 770 7 at low shear 
rates to a downward curved piece, followed by a point 
of inflection, and an upward curved piece, which finally 
goes over into a second linear behavior at very large shear 
rates, where ct = ry^^ 7. This S-shape can be recognized 
in Figs. [3] and m Because of the finite slope of log^g a ver- 
sus logj^Q 7 at the point of infiection, one may speculate 
about an effective power-law logj^o o ~ c + c' logj^g 7- 
Fig. [3] this happens at Peo ~ 10^^. Yet, the power-law 
is only apparent because the point of infiection moves, 
the slope changes with distance to the glass transition, 
and the linear bit in the fiow curve never extends over an 
appreciable window in 7;^ 

Non-trivial power-laws in the fiow curves exist close to 
the transition itself. At e = 0, a generalized Hershel- 
Bulklcy law holds 

a(7 - 0) ^ at f 1 + 1^1" + C2 l^^p™ + cg l^l^") 
V 7* 7* 7* / 

(30) 

and describes the fiow curves over an appreciable part of 
the range 7 < 7* , where structural relaxation dominates 
the stress j^^'^^ the exponent is to = 0.143 in the F^2''" 
model for this A. It provides a semi-quantitative fit of 
the flow curves for more than a decade in 7 close to the 
glass transition as shown in Fig. [51 There |e| is quite small 
at these effective packing fractions. A qualitative differ- 
ence of the glass flow curves to the fluid S-shaped ones, 
is that the shape of cr{'^) constantly has an upward cur- 
vature in double-logarithmic representation. The yield 
stress can be read off by extrapolating the flow curve 
to vanishing shear rate. In Fig. [5] this leads to a value 
cr+ « 0.24 kBT/R% at (j)^^ = 0.622, which is in agree- 
ment with previous measurements in this system over 
a much reduced window of shear rates.— While this 
agreement supports the prediction of an dynamic yield 
stress in the ITT approach, and demonstrates the useful- 
ness of this concept, small deviations in the flow curve at 
low 7 are present in Fig. [5l We postpone to Sect. VI. E 
the discussion of these deviations, which indicate the ex- 
istence of an additional slow dissipative process in the 
glass. Its signature is seen most prominently in the loss 
modulus G"{uj) in Fig. El 

The storage modulus of the glass shows striking elastic 
behavior. G'{uj) exhibits a near plateau over more than 
three decades in frequency, which allows to read off the 
elastic constant Goo easily. 



C. Microscopic MCT results 

Included in figures [31 to [51 arc calculations using the mi- 
croscopic MCT given by Eqs. pTT)) to evaluated for 
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FIG. 6: Flow curves for 0efl = 0.580 (diamonds) and 
4>c?i = 0.608 (circles) rescaled with the parameters from the 
respective F^j^'^ioclel fits; the corresponding fits are given 
as red solid lines, with separation parameter e = —2 10~* 
((jicff = 0.580) and e = 310"" (ijiefi = 0.608). The blue 
dashed cuves gives the generalized Hershel-Bulkley law from 
Eq. (|30|l . which holds at e = 0; for the Fj^'-model, its param- 
eters are 7* = 5.06 lO^'T = 2.16 10^''/fo (corresponding to a 
Peclet number PeJ = 4.25 10"=* at 0eff = 0.608), ca ~ 0.936, 
and C3 « 0.632. The green solid curve is the (critical) fiow 
curve at £ = . The magenta dashed straight line labeled rfl^ 
denotes the high-shear asymptote Eq. (|27p . 

hard spheres 1^ This is presently possible without shear 
only (7 = 0), because of the complications arising from 
anisotropy and time dependence in Eq. (|12p . The a pri- 
ori unknown, adjustable parameter is the matching time 
scale which we adjusted by varying the short time 
diffusion coefficient appearing in the initial decay rate 
in Eq. ([Tl]) . The computations were performed with 
rq(t) = Fq = Dgq^/Sq, and values for Ds/Dq are re- 
ported in the captions of Figs.[31to[ni and in table I. The 
viscous contribution to the stress is again mimicked by 
including 7700 like in Eq. (|25p . 

Gratifyingly, the stress values computed from the mi- 
croscopic approach arc close to the measured ones; they 
are too small by 40% only, which may arise from the 
approximate structure factors entering the MCT calcu- 
lation; the Percus-Yevick approximation was used herci^ 
In order to compare the shapes of the moduli the MCT 
calculations were scaled up by a factor Cy = 1.4 in Figs.[31 
to\5\ Microscopic MCT also does not hit the correct value 
for the glass transition point^i^ It finds cj)^'-^'^ = 0.516, 
while our experiments give 4'°^'^ w 0.58. Thus, when 
comparing, the relative separation from the respective 
transition point needs to be adjusted as, obviously, the 
spectra depend sensitively on the distance to the glass 
transition; the fitted values of the separation parameter 
e are included in Fig[51 

Considering the low frequency spectra in G'{uj) and 
G"{u!), microscopic MCT and schematic model provide 
completely equivalent descriptions of the measured data. 
Differences in the fits in Figs. [31 to [51 for uuIijj/Do < 1 
only remain because of slightly different choices of the 
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FIG. 7: The plots demonstrate that the reduced flow curves 
and the reduced moduli are unique functions only depending 
on (jtcff- AH flow curves are down curves. The flts using the 
schematic F^j'-model were performed with the data points 
at 13.01wt% taken before the onset of crystallisation (data to 
the right of the vertical bars) . Black diamonds: 12.10wt% and 
(/)oft = 0.527. Black circles: 13.01wt% and 0cff = 0.527. Red 
diamonds: 12. 10wt% and -^eff = 0.578. Red circles: 13.01wt% 
and (pad ~ 0.580. Green diamonds: 13.01wt% and (pcH ~ 
0.608. Green circles: 13.58wt% and (fj^g = 0.606. 



fit parameters which were not tuned to be close. These 
differences serve to provide some estimate of uncertain- 
ties in the fitting procedures. Main conclusion of the 
comparisons is the agreement of the moduli from micro- 
scopic MCT, schematic ITT model, and from the mea- 
surements. This observation strongly supports the uni- 
versality of the glass transition scenario which is a central 
line of reasoning in the ITT approach to the non-linear 
rheology. 

At large 7 and large w hydrodynamic interactions be- 
come important. In the flow curves, ri2^, and, in the 
loss modulus, rj'^ become relevant parameters, and the 
structural relaxation captured in ITT and MCT is not 
sufficient alone to describe the rheology. Qualitative dif- 
ferences appear in the moduli, especially in G"(tt'), be- 
tween the schematic model and the microscopic MCT. 
While the storage modulus of the F^^^-model crosses over 
to a high-w plateau already at rather low w, the mi- 
croscopic modulus continues to increase for increasing 



frequency, especially at lower densities; see the region 
Lu ^ 10'^ Do / Rfj in Figs. [3] to m The latter aspect is con- 
nected to the high-frequency divergence of the shear mod- 
ulus of particles with hard sphere potential^ as captured 
within the MCT approximation /n.ft,^.^, As carefully dis- 
cussed by Lionbcrgcr and Russel, lubrication forces may 
suppress this divergence and its observation thus depends 
on the surface properties of the colloidal particlesi^ 
Clearly, the region of (rather) universal properties aris- 
ing from the non-equilibrium transition between shear- 
thinning fluid and yielding glass is left here, and particle 
specific effects become important. 

D. Parameters 

In the microscopic ITT approach from Sect. II, the 
rheology is determined from the equilibrium structure 
factor Sq alone. This holds at low enough frequencies 
and shear rates, and excludes the time scale parameter to 
of Eq. p?|) , which needs to be found by matching to the 
short time dynamics. This prediction has as consequence 
that the flow curves and moduli should be a function 
only of the thermodynamic parameters characterizing the 
present system, viz. its structure factor. 

Figure [7] supports this claim by proving that the rhe- 
ological properties of the dispersion only depend on the 
effective packing fraction, if particle size is taken account 
of properly. Figure [7] collects flow curves and moduli 
measured for different concentrations of particles accord- 
ing to weight, and for different radii Rh adjusted by 
temperature. Whenever the effective packing fraction, 
4>cS ~ (47r/3)n_R|j, is close, the rheological data overlap 
in the window of structural dynamics. Obviously, ap- 
propriate scales for frequency, shear rate and stress mag- 
nitudes need to be chosen to observe this. The depen- 
dence of the vertices on Sq (Eqs. (|10ll3p ) suggests that 
ksT sets the energy scale as long as repulsive interactions 
dominate the local packing. The length scale is set by the 
average particle separation, which can be taken to scale 
with Rh in the present system. The time scale of the 
glassy rheology within ITT is given by to from Eq. (|17p . 
which we take to scale with the measured dilute diffusion 
coefficient Dq. Thus the rescaling of the rheological data 
can be done with measured parameters alone. Figure [7| 
shows quite satisfactory scaling. Whether the particles 
are truely hard spheres is not of central importance to 
the data collapse in Fig. [7] as long as the static structure 
factor agrees for the (pcS used. Fits with the F^^^-model 
to all data are possible, and arc of comparable quality to 
the fits shown in Figs. [3] to [5] 

The fitted parameters used in the schematic F^2''" 
model are summarized in Fig.[8l Parameters correspond- 
ing to identical concentrations by weight are marked by 
identical colours. Within the scatter of the data one may 
conclude that all fit parameters depend on the effective 
packing fraction only. This again supports the mentioned 
dependence of the glassy rheology on the equilibrium 
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FIG. 9: Newtonian viscosity 770 (diamonds, left axis), elas- 
tic constant Goo (squares), and yield stress (t+ (circles; data 
rescaled by a factor 30; both Goo and o-"*" right axis), as 
functions of the effective packing fraction <j)cs as obtained 
from the fits performed with the F^^' -model. Filled sym- 
bols indicate where direct measurements of r^o were possi- 
ble. Black symbols: 10. 75wt%, red symbols: 12. 10wt%, green 
symbols: 13.01wt%, blue symbols: 13.58 wt%. The line gives 
a power-law fit to the viscosity-data over the full range us- 
ing the known 7 = 2.34 exponent from MCT, logj^g ?7o = 
A — 7 ■ logjQ {(j>1s — (jtcft)', the critical packing fraction is found 
as (ji^ff — 0.580. Horizontal bars denote the critical elastic 
constant GJ^ , and the critical yield stress cr^ as extrapolated 
from Fig. |6l 



FIG. 8: The fitted parameters of the Fjj -model (open sym- 
bols). Black symbols: 10.75wt%, red symbols: 12.10wt%, 
green symbols: 13.01wt%, blue symbols: 13.58 wt%. e and 
7c are dimensionless. Filled magenta symbols, included in the 
upper inset, give the e values fitted in the microscopic MCT 
calculations for 13.01wt%. The unit of Vc is ksT/R^ while 
r is given in units of Do/Rjj. The high frequency and high 
shear viscosities rj'^'^ are given in units of ksT /{DoRh)- 



structure factor. The initial rate F, which sets t^, appears 
a unique function of 4>es , also; an observation which is not 
covered by the present ITT approach. It suggests that 
hydrodynamic interactions appear determined by 0off in 
the present system also. 

Importantly, all fit parameters exhibit smooth and 
monotonous drifts as function of the external thermo- 
dynamic control parameter, viz. c/ieff here. Nevertheless, 
the moduh at low frequencies (e.g. G"(w) at tuRjj/Do = 
0.01), or the stresses at low shear rates (e.g. ct(7) at 
•jRljj/Do = 10~^) change by more than an order in mag- 
nitude in Figs. [5] to [51 Even larger changes may be ob- 
tained from taking experimental data not shown, whose 
fit parameters are included in Fig. [S] It is this sensi- 
tive dependence of the rheology on small changes of the 
external control parameters that ITT addresses. 

When comparing the parameters from the schematic 
model to the ones obtained from the microscopic MCT 
calculation of the moduli, one observes qualitative and 
semi-quantitative agreement; see the captions to Figs. [3] 
to [51 table I, and the upper inset of Fig. [51 For example, 
the increase of the prefactor of stress fluctuations is 
captured in the microscopic vertex where Sq enters (this 



follows because the rescaling factor Cy is density inde- 
pendent). Also the hydrodynamic viscosity rjoo = r]'^ 
roughly agrees and may be taken (/)off -independent in the 
fits with the microscopic moduli. The fitted values of 
r]oc are actually not too diiferent from data obtained in 
Stokesian dynamics simulations of true hard spheres, sup- 
porting our simplified view on the particle interactionsi^ 
On closer inspection, one may notice that the separa- 
tion parameter of the microscopic hard sphere calcula- 
tion obtains larger positive values than e fitted with the 
schematic model. Moreover, it follows an almost linear 
dependence on the efiective packing fraction as asymp- 
totically predicted by MCT, e w 0.65 (0cff - 0off)/'^off 
with glass transition density (pc = 0.587 slightly higher 
than from the schematic model fits. The differing be- 
havior of the separation parameter from the fits with the 

F^2'*"niodel in the glass is not understood presently. The 
microscopic calculation signals glassy arrest more clearly 
than the schematic model fit. The short time diffusion co- 
efficient Da/ in the microscopic calculation decreases 
as expected from considerations of hydrodynamic inter- 
actions. Gratifyingly we find values in the range of the 
short-time self diffusion coefficient observed in Stokesian 
dynamics simulations for hard spheresiSi The initial rate 
r, however, of the schematic model increases with pack- 
ing fraction. The ad hoc interpretatation of F as mi- 
croscopic initial decay rate evaluated for some typical 
wavevector g*, viz. the ansatz F = Dsq^/Sq_^, thus ap- 
parently does not hold. 

While the model parameters adjusted in the fitting 
procedure only drift smoothly with density, the rheo- 
logical properties of the dispersion change dramatically. 
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TABLE I: Parameters of the fits with the microscopic MCT 
to the linear-response moduli G' [uj) and G"{lo). The first two 
columns of separation parameter e and short-time diffusion 
coefficient ratio Ds/Dq correspond to the fits shown in Figs.|3] 
to [5] and Fig. 1111 fsolid lines), while the second columns of e' 
and D'^/Dq correspond to the dashed- lines in Fig. 1111 when 
no value is given, the values from the first two columns apply. 
In all cases Cy = 1.4 and 7700 = 0.3 ksT / (DqRh) are used. 

Figure [9] shows the Newtonian viscosity as obtained from 

extrapolations of the fits in the ¥[2 -model. It changes by 
6 orders in magnitude. From the combination of G"{uj)- 
and flow curve data we can follow this divergence over 
more than one decade in direct measurement. From the 
divergence of r]Q the estimate of the critical packing frac- 
tion can be obtained using the power-law Eq. (HH]), be- 
cause the exponent 7 is known. We find (j)'^ff = 0.580 in 
nice agreement with the value expected for colloidal hard 
spheres. On the glass side, the elastic constant and yield 
stress jump discontinuously into existence. Reasonable 
values are obtained from the F^2''"niodel fits compared 
to data from comparable systems. The strong increase 
of the elastic quantities upon small increases of the den- 
sity is apparent. 



E. Additional dissipative process in glass 

One of the major predictions of the ITT approach con- 
cerns the existence of glass states, which exhibit an elastic 
response for low frequencies under quiescent conditions, 
and which flow only because of shear and exhibit a dy- 
namic yield stress under stationary shear. Figure [5] shows 
such glassy behavior, as is revealed by the analysis using 
the schematic and the microscopic model. Nevertheless, 
the loss modulus G"{lo) rises at low frequencies, clearly 
indicating the presence of a dissipative process. It is not 
accounted for by the present theory. Also, the storage 
modulus G'{lo) shows some downward bend at the low- 
est frequencies. 

These deviations can not be rationalized by ageing ef- 
fects or non-linearities in the response; see Fig. [101 We 
checked the dependence on time since quench to this glass 
state and also the linear dependence of the stress on the 
shear amplitude. While we find ageing effects shortly 
after cessation of pre-shearj^^i^i^i these saturate after 
one day, when the drifts of the spectra come to a stop. 
Ageing effects do not change the spectra qualitatively, as 
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FIG. 10: The storage G'{u) and loss G"(lj) moduU for dif- 
ferent waiting times The data have been also plotted as 
function of ut as suggested recently.— See text for further 
explanation. 



the dissipative process appears to possess a finite equi- 
librium relaxation time. As suggested recently for dense 
PNIPAM microgel dispersions^^ the same data have been 
plotted as function of uot. Here, t is the total waiting 
time and is defined as function of the waiting time tw 
before starting the measurement and the time 5{t{u)n)) 
expired between and the acquisition of the data as 
t = t{iOn) =tw + S{t{u}n))- The curves collapse on a mas- 
ter curve in the low frequency range up to uut w 3000 as 
expected from ageing theory for waiting time < 8200 s. 
This prediction is no more respected for longer waiting 
times, where an additional relaxation process is identi- 
fied. It cuts off the ageing behavior, when the age of the 
sample approaches the value of its relaxation time. This 
supports the introduction of an hopping phenomenon in 
our model, with a characteristic relaxation time of the 
order of lO'' s Rj^/SDo = lO^Rjf/TDo = 8.8 10^ s see 
FigO. 

Let us stress, moreover, that the state shown in Fig. [5] 
is not a fluid state within the present approach. The 
presence of an elastic window in G'(w), its increase as 
function of packing fraction, and the upward curvature 
of the flow curves rule out a negative separation pa- 
rameter e < of this state at (/)off = 0.622. Calcula- 
tions within the microscopic MCT document this con- 
vincingly. Figure [TT] compares the MCT calculations for 
hard spheres with moduli ranging from fluid to glassy 
states. By adjusting the effective packing fraction, MCT 
semi-quantitatively describes the dominating modulus, 
either loss or storage one, for all states (corresponding 
curves already shown in Figs [3] to [5]) . At high concentra- 
tions, it describes the storage modulus G'{lo) on an error 
level of IksT / Rfj, and misses the loss modulus G"{lo) by 
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FIG. 11: Fits with microscopic MCT to the linear-response 
moduli G'{ui) (upper panel) and G" (uo) (lower panel) for the 
packing fractions (j)aa ~ 0.527 (black diamonds and lines), 
0eff = 0.540 (violet), (^cs = 0.567 (light blue), (p^s = 0.580 
(red), (j>cS ~ 0.608 (dark blue), and 4>cs — 0.622 (orange). 
Continuous lines give the fits optimized for describing the 
storage modulus G'{ijj); these fits are also shown in Figs. |3] to 
O and the corresponding parameters are included in Fig. [S] 
and summarized in the left two columns in Table 1. Broken 
lines for cpaH = 0.580 (red, overlapping with the solid light blue 
curve), and (fics = 0.608 (dark blue; the same curve would fit 
4>cff ~ 0.622) show microscopic MCT calculations attempting 
to fit the minima in G" {u) enforcing negative separation pa- 
rameters e (parameters included in table I). These fluid like 
spectra can rationalize G"{lo), but fail qualitatively to de- 
scribe G'{uj). Vertical bars in corresponding colours denote 
the frequencies below which crystallisation affects the data at 
the different (pcft- 



a similar absolute error. Yet, because the latter is itself 
of the order of G" (uj) « IkBT/R^ in the measurements, 
MCT fails to describe G"(tj) adequately. If, however, 
the effective packing fraction in the MCT calculations 
is adjusted to match the loss modulus G"{uj), then this 
fit fails completely to capture G'{uj) at high densities; see 
the dashed lines in Fig.[TT] Because the storage modulus, 
however, dominates the linear mechanical response of the 
glass, the second fit needs to be rejected. In conclusion, 
MCT correctly identifies the transition to a glass at high 
densities with dominating elastic response and yielding 
behavior under flow. It misses an additional dissipativc 
process, which contributes on the 10% level to the shear 
moduli and stresses in the frequency and shear rate win- 
dow explored in our experiments. 

The existence of an additional dissipative process con- 
tradicts the notion of 'ideal' glass states as described by 
the present ITT or MCT aproach. Clearly, the system 
at (/)efi = 0.622 becomes a fluid at even longer times. 



or lower shear rates and frequencies than observed in 
Fig. O This does not, however, contradict the obser- 
vation that the structural relaxation as captured in the 
ITT equations has arrested. In an extension of the ITT 
approach, it is possible to account for the additional de- 
cay channel in a an extended schematic model; sec the 
Appendix for details. Results from this extended F^2''" 
model arc included in Fig. [5] and demonstrate that none 
of the qualitative features discussed within ITT change 
at finite frequencies or shear rates. The additional pro- 
cess leads to fluid behavior at even lower uj or 7, and 
needs to be taken into account only, if exceedingly small 
frequencies or shear rates are tested; its relaxation time 
at (/)off = 0.622 exceeds lO^Rj^/Do 8.8.10^ s. It does 
not shift the location of the 'glass transition' as defined 
within the idealized ITT (MCT) approach, because this 
is already determined by the shapes of the flow curves 
and spectra in the observed windows. 



VII. CONCLUSIONS 

In the present study, we explored the connection be- 
tween the physics of the glass transition and the rheology 
of dense colloidal dispersions, including in strong shear 
flow. Using model colloidal particles made of thermoscn- 
sitive core-shell particles, wc could investigate in detail 
the vicinity of the transition between a (shear-thinning) 
fluid and a (shear-molten) glass. The high sensitivity of 
the particle radius to temperature enabled us to closely 
vary the effective packing fraction around the critical 
value. We combined measurements of the equilibrium 
stress fluctuations, viz. linear storage and loss moduli, 
with measurements of flow curves, viz. nonlinear steady 
state shear stress versus shear rate, for identical exter- 
nal control parameters. In this way we could verify the 
consequences from the recent suggestion, that the glassy 
structural relaxation can be driven by shearing and in 
turn itself dominates the low shear or low frequency rhe- 
ology. 

In the employed theoretical approach, the equilibrium 
structure as captured in the equilibrium structure fac- 
tor Sq sufficed to describe all phenomena qualitatively. 
As only exception, we observed an ultra-slow decay of 
all glassy states that is yet not accounted for by theory. 
Microscopic calculations were possible for the linear re- 
sponse quantities using mode coupling theory applied to 
hard spheres. Schematic model calculations were possible 
within the integration through transients approach, and 
simultaneously captured the linear and nonlinear rheol- 
ogy using identical parameter sets. Semi-quantitative 
agreement between microscopic and schematic model cal- 
culations and with the measured data for varying effec- 
tive packing fraction could be achieved adjusting a small 
number of fit parameters in smooth variations. 
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APPENDIX A: EXTENDED MODEL 
INCLUDING HOPPING 

The ITT equations contain the feed back mechanism 
that the friction increases because of slow structural rear- 
rangements. In the schematic F^^^-model this is captured 
by the approximation for the generalized friction kernel 
m{t) in Eq. ((22)) . For 7 = it leads to non-ergodic glass 
states at large enough vertices vi^2- A dissipative process 
explaining the fiuidity of glassy states should renormal- 
ize the diffusion kernel A(t). Moreover, this mechanism 
should become more important the longer the relaxation 
time in m{t). If, however, the additional dissipative pro- 
cess is too strong, all effects of the bare ITT approach 
are smeared out and the described phenomenology of the 
glass transition can not be observed. 

Gotze and Sjogren found when considering (possibly 
unrelated) dissipative processes in simple liquids that this 
can by achieved by splitting the diffusion kernel into two 



decay channels, one connected to the original m(t), and 
the other one connected to the new dissipation mecha- 
nism. In order for the second decay channel to take over 
in glassy states, it suffices to model it by one additional 
parameter (5 in a linear ansatz A'^'*''^'P(i) = 5m{t). This 
leads to the following replacement of Eq. (PT|) in the 

F<i^ -model 

at$(i) + r|$(t) +jdt'm{t-t') = o 

(Al) 

The memory function m{t) is still given by Eq. (|22[) 
because shearing decorrclates arbitrary fluctuations via 
shear-advection. All parameters of the model are kept 
as specified in Sect. IV, and solutions of this extended 
model with parameter 5 given in the caption arc included 
in Fig. [5] Importantly, the fluid like behavior in the rlie- 
ology at exceedingly small 7 and oj can now be captured 
without destroying the agreement with the original ITT 
at higher parameters. Apparently, a single parameter 5 
is not sufficient to model the non-exponential shape of 
the final relaxation process in the glass. Yet, further 
extensions of the model in order to describe this non- 
exponentiality go beyond our present aim. 
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